import xarray as xr
ds = xr.open_dataset('era5_data/surface_2023-01-01.nc')
print("surface的所有变量：",ds.data_vars)  # 打印所有变量
# 运行目录需要放在根目录
print("变量最大值统计：")
for var_name in ds.data_vars:
    var = ds[var_name]  # 取出变量
    max_val = var.max().item()  # 最大值（转为 Python 标量）
    shape = var.shape
    print(f"{var_name:<10} shape={shape}  max={max_val:.4f}")

# ds = xr.open_dataset('era5_data/upper_2023-01-01.nc')
# print("upper的所有变量：",ds.data_vars)  # 打印所有变量
